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We measure the spectrum of low- lying nucleon resonances using Bayesian fitting methods. We compare the 
masses obtained in the quenched approximation to those obtained with two flavours of dynamical fermions at 
a matched lattice spacing. At the pion masses employed in our simulations, we find that the mass of the first 
positive-parity nucleon excitation is always greater than that of the parity partner of the nucleon. 



The observed spectrum of excited nucleon 
states contains some puzzles, such as the anoma- 
lously light masses of the Roper AT* (1440) and 
A(1405). Neither of these states naturally lies 
within the quark model, and this had led to spec- 
ulation as to whether they are really three quark 
states. This spectrum, at least in principle, is 
accessible to lattice QCD calculations. In this 
study we extract the masses of the nucleon (N + ), 
its parity partner (N~) and the lightest positive- 
parity excitation (A' + ) and examine their order- 
ing. 

Lattice determinations of the nucleon mass 
generally use the following three interpolation op- 
erators: 
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These operators have an overlap with both 
positive- and negative-parity states, but on a lat- 
tice with (anti-)periodic boundary conditions we 
can use the parity projection operator |(1 ± 74) 
to delineate forward and backward propagating 
states of opposite parities. The time dependence 
of the parity projected correlators is then 



c±(t)=Y,Ai 



-Aff=t 



(4) 



The "diquark" part of N\ ; 3 couples the upper 
spinor components, whilst the N2 operator cou- 
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Figure 1. Mass versus number of exponentials. 
p = 6.2, k = 0.1346. 

pies upper and lower components, and so van- 
ishes in the non-relativistic limit. The expecta- 
tion, confirmed by the data, is that the A1.3 oper- 
ators give a much better overlap with the nucleon 
ground state while the N 2 operator has an overlap 
with the lowest positive-parity nucleon excitation. 

For the negative-parity states the signal is short 
lived Q, so the choice of fit range may lead to 
a substantial systematic uncertainty. The aim 
of Bayesian fitting is to eliminate this uncer- 
tainty Qj. Furthermore, by fitting all the data, it 
may be possible to measure the contribution from 
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Figure 2. Quark mass dependence of baryon 
masses at (3 = 6.2 

excited states. The number of exponentials A cxp 
is chosen such that the results for the ground and 
first excited states are stable if A exp is increased. 
The ground-state priors were guessed from the 
correlator data at large temporal separations us- 
ing the effective mass, 

(C{t + 1) 



M ett {t) = In 



(5) 



V c(t) 

and priors for the mass splitting and excited 
state amplitudes obtained using AM = Mf(r)- 
M{j s (t), where r < t, and Mf ff is determined from 
the correlator C with the ground state prior con- 
tribution subtracted. This ensures that the min- 
imisation algorithm starts in the neighbourhood 
of the global minimum of \ 2 . The width of the 
priors is chosen to be larger than the splitting. 
Furthermore, the fits are more stable if both M 
and AM are constrained to be positive, so instead 
of fitting to eqn we fit the following expression 
to the data 



f(t) =^A„exp 
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where X k = ln/x fe , fi a = m and fi k = M k -M k ^ X - 
The calculation was performed using the 
standard Wilson gluon action, and the non- 
perturbatively C(a)-improved Sheikholeslami- 
Wohlert action. Lattice artefacts thus appear in 



Table 1 

Details of the data sets, a -1 in GcV, the m^/rrip 
values for the DF configurations are at ft sca = 



,val 



cfgs 



V 



6.2,0 
6.0,0 
5.93,0 



216 2.91 
496 2.12 
623 1.90 



24 3 x 48 
16 3 x 48 
16 3 x 32 



N/A 
N/A 
N/A 



5.2,0.1350 
5.26,0.1345 
5.29,0.1340 



202 
102 
101 



1.91 
1.90 
1.94 



16 3 x 32 0.70(1) 
16 3 x 32 0.78(1) 
16 3 x 32 0.84(1) 



the spectrum at 0(a 2 ). The values of {3 and « sca 
for the coarsest quenched (Q) and the dynamical- 
fermion (DF) ensembles were chosen so that the 
lattice spacing is the same. The light hadron 
spectrum has already been determined by the 
UKQCD collaboration for [3 = 6.2 and (3 = 6.0 g 
and for DF matched ensembles j|. The N~ 
masses have also been determined for (3 = 6.2 
and $ = 6.0 §. 

We perform simultaneous fits to both the local- 
local and smeared-local correlators of both pari- 
ties. However, for the lightest two DF ensembles 
the limited number of configurations only allows a 
fit to the local-local correlator. Thus we only fit to 
one correlator for all the DF ensembles. The de- 
pendence of the extracted masses on the number 
of exponentials A^xp is shown in figure [l|, together 
with the priors, shown as the open symbols. The 
results are stable when A oxp > 10, corresponding 
to five states of each parity. For (3 = {6.2, 6.0} 



we use N r 



12, and for the DF and matched 



Q ensembles we use N e: 



10. 



The spectrum as determined from the opera- 
tors Ai } 3 (Ni hereafter) and N2 is shown in fig- 
ure |^; we find consistency between the masses 
of N~ using the two operators, and between the 
excited-state mass N[ + and ground-state mass 
A2". However, as the quark mass is decreased 
the data quickly become very noisy, and at lighter 
quark masses, the masses from the N2 operator 
are larger than those from the Ai. Only a sim- 
ple linear chiral extrapolation is attempted as the 
data is too noisy for a more sophisticated ap- 
proach. 

The comparison between the Q and DF ensem- 
bles at matched lattice spacing is shown in figure 
|^, revealing no systematic difference between the 
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Figure 3. Nf = 2 (open symbols) versus Nf = 
(closed symbols) . Some of the symbols have been 
offset for clarity. 

Q and DF results. In all cases the mass of N' + 
is above that of N~. Furthermore the positive- 
parity excited state mass obtained from N\ is con- 
sistent with the ground-state mass obtained from 
N 2 , except at light quark masses where the lat- 
ter becomes very noisy. Finally, the lattice spac- 
ing dependence of the measured masses on the 
quenched ensemble is shown in figure ||. 

We have measured the nucleon spectrum and 
found that the masses are ordered according to 
N'+ > N~ > N+ for both the Q and DF data 
at all measured quark masses and lattice spac- 
ings, thus providing little evidence for the Roper 
resonance as a simple three-quark nucleon exci- 
tation. However, there are several caveats. The 
statistical quality of the data is poor. The spatial 
extent is only around 1.7 fm, though a compar- 
ison with a calculation of both the N + and N~ 
masses on a lattice of around twice that volume 
at (3 = 6.0 suggests that at these quark masses 
the finite- volume corrections are small |l| ; such a 
comparison has not been performed for N' + . Fi- 
nally, the quark masses are still large, with the 
smallest DF mpsjmy — 0.7. A calculation in 
the quenched approximation to QCD using over- 
lap fermions with pseudoscalar masses as light as 
200 McV indicates a dramatic cross-over between 
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Figure 4. Lattice spacing dependence in quenched 
approximation to QCD. 

N~ and N /+ at light quark masses jq]. 

The experimental study of the nucleon reso- 
nances is a crucial tool for investigating the na- 
ture of the QCD interaction. Lattice calculations 
have a vital role both in interpreting the data, 
and in guiding future experimental searches. 
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